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Abstract 

The linearly polarized gravitational waveforms from a certain type of ro- 
tating, evaporating cosmic string - the Allen-Casper-Ottewill loop - are 
constructed and plotted over the lifetime of the loop. The formulas for 
the waveforms are simple and exact, and describe waves which attenuate 
self-similarly, with the amplitude and period of the waves falling off linearly 
with time. 

Short Title: Waveforms from a String Loop 
PACS numbers: 04.25.Nx, 98.80.Cq 

1. Introduction 

Cosmic strings are thin filaments of topologically-trapped Higgs field energy which 
may have formed at a symmetry-breaking phase transition in the early Universe 
(see PP for a review). In the zero-thickness approximation they are effectively 
line singularities whose dynamics and stress-energy content are governed by the 
Nambu-Goto action. Knowledge of the stress-energy tensor can in turn be used to 
model the string's gravitational field, and in particular the effects of gravitational 
back-reaction on the string's own motion. 



In two previous papers [21 El I examined the back-reaction problem for a cer- 
tain type of cosmic string loop, the AUen-Casper-Ottewill (AGO) solution [1], 
which in the absence of gravity is rigidly rotating, and is thought to have the 
lowest gravitational radiative efficiency - and therefore the longest lifetime - of 
any flat-space loop configuration. In these papers, I constructed a solution of the 
linearized Einstein equations describing the self-similar evaporation of the AGO 
loop, which radiates energy and angular momentum at a constant rate and evap- 
orates completely in a finite time. The solution can also be matched to a relic 
vacuum spacetime across the future light cone of the final evaporation point of 
the loop. 

Given the complex nature of the Einstein equations for an extended object like 
a cosmic string, it is unlikely that a fully non-linear solution to the back-reaction 
problem will ever be developed in even the simplest cases, but fortunately it 
can be shown that the self-gravity of a GUT string would almost everywhere be 
small enough to justify a weak- field treatment [5]. There are thus good reasons for 
believing that the evaporating AGO loop described in [5] very closely approximates 
the evolution of a realistic loop configuration, if indeed cosmic strings ever existed 
in the early Universe. 

In the present paper I continue the analysis of the evaporating AGO loop 
by calculating and plotting the linearly polarized waveforms emitted by this loop. 
Since the waveforms emitted by an object are discernible only in the object's wave 
zone (at distances large compared to characteristic scales within the object), all 
the information needed to construct them can be extracted from the weak-field 
limit. The resulting analytic expressions for the waveforms from the evaporating 
AGO loop - which are derived in detail below - are unexpectedly simple, and are 
evidence once again of the tractability of the solution. 

Allen and Ottewill in an earlier study of gravitational wave spectra from string 
loops [6] have calculated and plotted the waveforms emitted by the stationary 
AGO loop, but the expressions they derived are - unlike the formulas developed 
below - infinite Fourier sums, and of course these omit the amplitude attenuation 
and frequency acceleration that accompany the evaporation. Although current 
limits on the mass per unit of cosmic strings make it very unlikely that spectra 
from individual loops will ever be observed, the calculations that follow offer an 
unusual example of a complete analytic set of waveforms from a realistic, extended, 
self-gravitating system. 
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2. The Evaporating AGO Loop 



The world sheet of a zero-thickness cosmic string is the two-dimensional surface it 
traces out as it moves, and is described parametrically by a set of equations of the 
form x°- = X"-{(^), where = [x^, x^, x^, x^] = [t, x, y, z] are local coordinates on 
the four- dimensional background spacetime (JsA^Qab)-, and the parameters C,^ = 
(C°,C^) are gauge coordinates. In terms of the intrinsic two-metric 

IaB = gabX^-^A X^,B (2.1) 

(with X'^,A shorthand for dX"'/d(^) the Nambu-Goto action [3 [8] reads 

/ = -/^ / l'^'d\, (2.2) 

where 7 denotes | det(7^^)| and /i is the (constant) mass per unit length of the 
string. 

The corresponding stress-energy tensor T"^ can be constructed from the first 
variation of / with respect to Qab in the standard way, and has the form 

T''\x') = fig-^/^ J (5^') (x" - X") (2.3) 

where g = \ det{gab)\ and p"* = ■y^^X"-,AX^,B is the projection tensor onto the 
tangent space of the world sheet. The equation of motion of the string is similarly 
constructed by setting the first variation of / with respect to X" to zero, and 
reads 

g^(7^^X^cz. +P"^"rtj = 0, (2.4) 

where g"^ = g^-^—p"-^ is the orthogonal complement of p"'^, and F^^ is the Christoffel 
symbol associated with gab- 

If the background spacetime is flat, so that gab = Vab where ?7„j, is the Minkowski 
metric tensor (here taken to have its rigid form 77 = diag (1,-1,— 1,-1) ), then 
the equation of motion (12. 4p reduces to the wave equation 

X°',rT = X"',a-a- (2.5) 

provided that the coordinates (r, a) = (C°, C^) are chosen to satisfy the gauge 
conditions X^ + X^ = and X^ ■ X„ = (which is equivalent to choosing 7^^ to 
be diagonal and trace-free). 
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A particular solution of the flat-space equation of motion fl2.5p is the AUen- 
Casper-Ottewill (AGO) loop which is described by the trajectory 



X' 



a 



J_ 



Lq[u + f , cosf , sinf , \u 




(2.6) 



where u = 27r(r + (t)/Lq and v = 2tt{t — <7)/Lq, the constant Lq is the invariant 
length of the loop, and the parameter r has been identified with the coordinate 
time t. The gauge coordinate u covers the range (— tt, tt], with the limit u — tt"*" 
coinciding with u = n, while the range of v is unrestricted. The AGO loop rotates 
rigidly with fundamental oscillation period tp = Lq/2, as shown in Figure 1, and 
has total energy E = fiLo and angular momentum J = -^iiLqZ. The points at 
the extreme top and bottom of the loop, which correspond to m = vr and u = 
respectively, are known as kinks. Here the tangent vector is discontinuous, 
although the trajectory itself remains continuous. 

Given any flat-space solution of the equation of motion (12. 5p . it is possible to 
calculate the gravitational power radiated by the loop by taking the Fourier trans- 
form T"'' of the stress-energy tensor (12. 3p and substituting it into the quadrupole 
formula for the power per unit solid angle at future null infinity: 



where = 27r/tp is the circular frequency of the source and m is the Fourier wave 
number. Because T"-^ is proportional to the string's mass per unit length n, the 
total power P of the loop scales as fi^, and a dimensionless measure of the power 
is the loop's radiative efficiency 7° = P/ /i^. The significance of the AGO loop lies 
in the fact that its radiative efficiency 7° ^ 39.0025 - which was first calculated 
by Allen, Gasper and Ottewill [4J in 1994 - is the lowest of any known cosmic 
string loop. All other known analytic and numerically-generated solutions have 
higher efficiencies [H O HO] , although a definitive proof that the AGO loop has 
the lowest possible value of 7^^ has not yet been found. 

The stress-energy tensor (12. 3p can also be used to estimate the gravitational 
field induced by any fiat-space loop by inserting T"'' into the linearized Einstein 
equation. In the linearized approximation, the metric tensor is decomposed in 
the form gat = Tj^f^ + hab, where the components of hah are assumed to be small. 
If hab is constrained to satisfy the harmonic gauge conditions h^^ih = |^,a (with 
h = hi) then at first order in hab the Einstein equation G"* = — SvrT"* reads 
Ohab = —IGTcSab, where Sab = Tab — ^Vab^c aiid O = of — is the fiat-space 




(2.7) 



m=l 
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d'Alembertian. The corresponding retarded solution for hat is: 




(2.8) 



with x"" = [t, x]" and t' = t — |x — x'| the retarded time at the source point x'. In 
view of the equation fl2.3p for T'^'', the components of hat are all of order /i. 

Once hab has been calculated, it is in principle possible to determine the grav- 
itational self-force on the loop, and then use the linearized version of the equation 
of motion (12 ■4p to model the effect of the back-reaction averaged over a single 
oscillation period of the loop. This has been done for the AGO loop in where 
it was found that after a single period At = tp the loop retains its original shape, 
but its length Lq shrinks by an amount AL = —^•y'^fiLQ ^ — 19.501yuLo while its 
rotational phase advances by an angle of about 38.92/i radians. In particular, the 
radiative efficiency 7" remains unchanged, and so the AGO loop radiates energy 
at a constant rate even in the presence of gravitational back-reaction. 

In [3| I used a simple geometric argument to demonstrate that if the AGO loop 
does evaporate self-similarly with a constant power P = 7*^//^ then its trajectory, 
to leading order in /z, will have the form 



where k = •y^/iAn) and = j^;^Lq is the total lifetime of a loop whose length 
at time t = is Lq. As in the case of the stationary loop (12.61) . the coordinate 
u covers the range (— 7r,7r]. The evaporating loop (12. 9p shrinks to a point at 
x° = \tL, 0]" in the limit as 00. Furthermore, in the neighborhood of any 
spacelike surface defined by t = ti{l — e~'^^^'), where is an arbitrary constant, 
the evaporating trajectory (12. 9p is isometric to the stationary trajectory (12. 6p to 
leading order in /i, save with Lq replaced by L^ = Lqc''^^'"* , and cosf and sint> 
replaced by cos(w* -|- v) and sin(w^, -|- v). The evaporating trajectory therefore 
satisfies the fiat-space equation of motion (12.51) to leading order in /x. 

As with the stationary AGO loop, it is possible to use the retarded integral 
(12. Sp to calculate the metric perturbations hab induced by the evaporating loop 
(12. 9p . The exact form of these perturbations will be presented in Section 3. The 
perturbations were used to develop the string equation of motion (12.41) to order 
in [3], where after a lengthy calculation it was shown that the trajectory of the 
evaporating AGO loop, correct now to order /i^, has the form 



X"(M,t;) = [ti,0]" + ^Loe-''^("+^)[. 



— -j^, cosf , sinf , 1^1 — iyr 



a 



(2.9) 



a 



(2.10) 
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with 

f{u,v) = e^^(")[cos(t;)x + sin(t;)y] + [{\u\ - ^n) + fxA{u)]z, (2.11) 

where the functions F, G (both periodic in u with period vr) and A (anti-periodic 
in u, with A{0) = —A{7i) = — 127r + Svr In 27r + 2 Si(27r)) are fixed up to a Poincare 
transformation and an arbitrary gauge transformation of A. In addition, the 
expected value of k = 7°/(47r) - namely 4[ln27r + 75. — Ci(27r)]/7r ^ 3. 103 7 (where 
7e is Euler's constant) - emerges as a consistency condition on the solution f l2.10p . 

It was further demonstrated in [3] that the metric perturbations hab can be 
matched smoothly (to order /i) to a Minkowski spacetime across the future light 
cone Fi = {t — ti = |x|} of the final evaporation point at = [t^, 0]°, and so the 
evolving AGO loop (12. 9p evaporates completely to leave behind a remnant vacuum 
spacetime, as would be expected. The gross features of the loop's evaporation are 
illustrated in Figure 2. 



3. The Far-Field Metric Perturbations 



If [t,x] is a general field point satisfying the constraint t — t^ < |x| (and so lying 
to the past of the null surface Fi), then the linearized metric perturbations hab 
due to the evaporating AGO loop (12. 9p have the self-similar form [3| 



hab{t,yi) = 4/i 



/+ + /_ + -C+-C- 1+ 
S+ + S^ + s+ 

/+ + /- -C4 



(3.1) 



ab 



where I±, S± and C± are functions of the rotation angle 

^(t, x) = -{Kfi)-^ ln[(|x| - t + tL)/tL\ 



(3.2) 



and the components of an image point = [t,x]'^ lying on the future light cone 
-^0 = = |x|} of the origin, and defined by 



t 



--\-t + tL)/tL] ^|x|, X = [{\:>c\ - t + tL)/tL] ^{xcosi> + ysinip) 

t + tL)/tL]^^{ycosip — xsmip) and z = — t + tL)/tL]~^z. 

(3.3) 
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To be specific, if [T, X, Y, Z] = At^Lq [t, x, y, z] and the functions Vq and Vi 



are defined implicitly by the equations 



T-Vo = [{X- cos Vof + {Y- sin Vof + {Z + hi 



and 



T - (tt + V^i) = [(X - cos Vxf + (F - sin V^f 
then /-t, S± and (7^ have the explicit forms 



[Z-l-K 



,21 1/2 



l2ll/2 



(3.4) 



(3.5) 



/+ = 


= ln[(x+ - - 


-27r)/(x+-K))], 




- = ln[(X- 


- ^o)/(X- 




, (3.6) 






- cos(x- 


^ + t/;) sin(x+ + V^) 




Si(x+ - 


- 27r) - 


Si(x+ 


-V^o) " 






sin(x+ 


+ V^) cos(x+ + V^) _ 




. Ci(x+ - 


Vx - 27r) - 


Ci(x+ 





and 



- cos(x_ + i)) sin(x„ + i)) 
sin(x_ + ^) cos(x_ + V') 



Si(x_ 
Ci(x_ 



K)) 
K)) 



Si(X- 
Ci(x- 



(3.7) 



(3.8) 

^ cos w dw 



where x± = T±{Z+^), and Si(x) = sinw dw and Ci(x) = — 

are the usual sine and cosine integrals. Note also that T = (X^ + Y'^ + Z^)^/^, as 
lies on the light cone Fq. 

The bars on 5-1- and C± have here been used to distinguish them from the cor- 
responding unbarred quantities introduced in [3j, which were defined analogously 
to (13. 7p and (13.81) except with ip = 0. Apart from this explicit '0 dependence, the 
metric perturbations defined by (13.11) are determined completely by the compo- 
nents of the image point x° alone. It is for this reason that the weak-field metric 
was described in [3] as exhibiting "rotating self-similarity" Q 



-'^In fact, it can be shown that the vector field 

fc" = [0, -y, a:, 0]'' - - iL, X, z]^ 

which is constructed by taking fc° = dx"" /dip (with the components of the image point a;" kept 
fixed), is a conformal KiUing vector field of the linearized metric gat = Vab + ^ab, as it satisfies 
the equation 

to linear order in fi. 
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The main purpose of this paper is to calculate the gravitational waveforms 
generated by the metric perturbations hab as measured at a field point at a fixed 
distance r = |x| from x = 0. It will henceforth be assumed that the field point 
lies in the wave zone, where r ^ Lq. To calculate the waveforms it is necessary 
to know the metric components only to linear order in /i, and these only to order 

It should be noted, however, that the quantity r — t appearing in the defi- 
nition (13.21) of and subsequent equations is not necessarily large. In fact, the 
combination {t — r)/tL is a convenient measure of observer time at the field point, 
and ranges from on the initial null surface Fq to 1 in the limit as the field point 
approaches the boundary null surface F^. 

In view of fl3.3p . therefore, the components of the image point are all of 
order r in the wave zone, and the equations (13.41) and (13. 5p defining Vq and Vi 
reduce to 

= r-^[x cos(\/o + ^)+y sm{Vo + ^) - f 2] (3.9) 

and 

Vi + n = r-^[xcos{Vi+^) + ysm{Vi+^) + fz] (3.10) 

respectively. Although the functions Vq and Vi typically do not converge as the 
field point approaches F^ (where ip ~^ 00), it is clear that |Vo| and \Vi + n\ are 
both bounded above by (1 + |7r^)^/^. 

As the asymptotic behavior of the metric perturbations hat is dominated by 
the quantities x+ and X-, the cases in which the field point does or does not lie 
on the 2;-axis (where x = y = and so T = \Z\) need to be considered separately. 

If the field point does lie on the z-axis, then (13. 9p and (I3.10p immediately 
reduce to Vq = —^C, and Vi = — tt, where ^ =sgn(z), while x± = ±{Z + 
However, two of the metric functions S± and C± are then undefined, which is 
simply an indication that the asymptotic expansions (13. 9p and (I3.10p are too 
coarse. 

Fortunately, it is possible to solve equations (13. 4p and (13.50 explicitly in this 
case to give [3] 

Vo = ^Z-[l + iZ + f )2] V and = ^Z - n - [1 + {Z - f )2]V'. (3.11) 

The metric perturbation hf^ is then identically zero, while the remaining, non-zero 
components have the asymptotic forms 

htt = Kx = hyy = = -2^L/r + 0(r"^), (3.12) 
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^/ite = -hzx = {cos ilj)i^L/r + 0{r ^) 



(3.13) 



and 



Chty = -h,y = {smil;)iiL/r + 0(r-'), (3.14) 

where L is the parametric length of the string at the time of emission of the 
waveform: 

L{t)^Lo[ir-t + tL)/tL], (3.15) 

In fact, as will be seen in the next section, the corresponding waveforms carry no 
gravitational energy along the 2;-axis. 

In the more general case in which the field point lies off the 2;-axis, the functions 
I±, S± and C± can be represented by the leading-order expressions 

/+ = (1 + cos^)-i(Vo -Vi- 27r)L/(47rr), (3.16) 

7_ = (1 - cos^)-^(yo - Vi)L/(A7rr), (3.17) 

l + cos^)-iL/(47rr) (3.18) 















and 








' s_ ' 










where cos 9 - 





cos(Vl + '?/')— cos(Vo + i^) 
- sm{Vi + V') + sm{Vo + ^j) 



- cos{Vi + cos(\/o + ij) 
sin(Vi + -0) — sin(\/o + "0) 



(1 - cos^)-^L/(47rr), (3.19) 



The components of the hab therefore have the asymptotic forms: 

2i2L 



htt = k 



XX hyy 



htz 
htx — 

hty = 



irr sin^ 9 
2/xL 



[-7r+ (Vi + TT-Vo) cos 9], 



Tir siv? 9 
2fj,L cos9 . 
Trr sm 9 
2fj,L cos 9 



{Vq -Vi-TT + 7rcos6'), 
cos{Vi + + cos{Vo + ip)], 



irr sin^ 9 



hzx -2/1 

Trrsm 9 



[- sm{Vi + tj;)+ sm{Vo + tp)], 
[cos{Vi + V') - cos(yo + ^)] 



and 



h 



zy 



2/xL 
7rr sin^ 9 



[sin(yi + -0) - sin(yo + V')]- 



(3.20) 

(3.21) 
(3.22) 
(3.23) 
(3.24) 

(3.25) 
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4. Calculating the Waveforms 



The procedure for calculating the linearly polarized gravitational waveforms 
and hx from a knowledge of hab has been conveniently summarized by Allen 
and Ottewill in [6J. This reference generates the waveforms for a number of 
simple cosmic string loop configurations, including the stationary AGO loop fl2.6p . 
However, the waveforms for the AGO loop are there presented in the form of series 
approximations, as the metric perturbations hab are calculated by means of Fourier 
expansions. By contrast, the waveforms constructed in this section are those from 
the evaporating AGO loop (12.91) . and these can be cast in a relatively simple closed 
form, as will be seen shortly. 

The first step in calculating the waveforms is to construct the radiative per- 
turbations hab formed by adding the potential term 2/iL/r to each of the diagonal 
elements of hab, so that 

hab = hab + i2fiL/r)6ab- (4.1) 

In the case where the field point lies on the z-axis this yields the very simple 
radiative metric 



h 



ab 





^ cosip 
^ sin?/' 




^ cosip 



— cosip 






— sinip 





cosip 
■ sinip 




(4.2) 



On the other hand, if the field point does not lie on the 2;-axis then 

2yuL cos 9 



htt = K 



h 



yy 



h. 



-(Vi + TT — Vo — TTCOSl 



(4.3) 



Tir sin^ 9 

whereas the remaining, off-diagonal components fl3.21l) - fl3.25l) are unchanged: 

hab = hab if a^b. (4.4) 

The next step is to eliminate the time-time and time-space components of hab 
by performing the gauge transformation 



h'ah ~ ^ab — ^aEb — nb£a, 



(4.5) 



where n° = [l,n]" is a null vector in the direction of the field point x", and 
Ea = [eo, e]a with 



hhtt 



(4.6) 
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and the components of the 3-vector e given by 

= -l^ttnj + htj. (4.7) 
If the field point is represented in standard spherical polar coordinates as 

= [t, r sin 6 cos 0, r sin 6 sin 0, r cos 6]^ (4-8) 

then 

n = (sin 6 cos (p, sin 6 sin 0, cos 6) (4-9) 

(and the quantities Uj are the components of — n). 

In the case where the field point lies on the z-axis, £0 = and n = (0,0,.^), 
while 

ej = {fiL /r){^ cos ij,^smilj,0)j. (4.10) 

It is easily verified that in this case h'^j^ = and so the metric perturbations are 
pure gauge terms. An immediate consequence is that the evaporating AGO loop 
radiates no gravitational energy along the z-axis, as mentioned earlier. 

In the more general case where the field point lies off the z-ax.is, it proves 
convenient to eliminate the term Vi + vr — Vq appearing in hu and ht^ by invoking 
equations (13.91) and (IS.lOp . which can be written more compactly as 



and 
Then 

and 

where 
Hence, 



Vo = cos(Vo + ^-0)sin^-f cos^ (4.11) 

+ TT = cos( + - 0) sin e + f cos ^. (4. 12) 

7 7 7 7 2/iLcOs6', , 

htt = h^x = hyy = = - 1^0) (4.13) 

htz = ^^(/ti - ^o), (4.14) 

vrr sm & 

Ko,i = cos{Vo,i+tp - (p). (4.15) 



^ -[{ki — kq) sin 6^ cos 6^, {{ki — kq) sin^ 9 cos0 — 2(ci — cq)} cos6', 



vrr sin^ 9 

{(ki — kq) sin^ 6^ sin — 2(si — Sq)} cos 9, — (/ti — kq){2 — cos^ 9) sin6']a, 

(4.16) 
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with 

co,i = cos(yo,i and so,i = sin(yo,i (4.17) 

The space-space components of the corresponding transformed metric pertur- 
bations h'^^ are: 

- ^ 2/iI.cosg ^^^^ _ ^^ ^^ ^ sin^^cos^^) - 2(ci - Cq) cos0], (4.18) 
vrr sm d 



- ^ 2/iLcos g - fi:o)(l + sin^ ^sin^ 0) - 2(si - sq) sin0], (4.19) 
'^'^ 7rr sm 



h' = — ^—{ki — Ko) sin6'cos^, (4.20) 
vrr 



= [-(/ti - /to) (1 + cos^ g) sin cos + (1 - 2 cos^ 0) (ai - (Tq)] , (4.21) 

^ vrr sm a 



and 



where now 



/i^^ = [(ci — Co) — (ki — Ko) sin^ 9 cos 0] (4.22) 
vrr 

Kz = -^[{si - So) - (ki - Ko)sin^6'sin0], (4.23) 



ao,i = sin(Vo,i+V^-0)- (4.24) 



An important feature of the 3-tensor /i^-^ is that the contraction n^h'ji^ is iden- 
tically zero, as can be verified directly from the expressions f l4.18l) - fl4.23l) . The 
metric perturbations h'jj^ can be reduced to canonical linearly polarized form by 
rotating the three spatial axes so that the new 2;-axis is aligned with the unit vec- 
tor n along the line of sight. This is accomplished by introducing an orthonormal 
triad {mi, m2, n} with 

m2 = (— sin0, COS0, 0) (4.25) 

and 

mi = m2 X n = (cos cos 0, cos ^ sin 0, — sin^^). (4.26) 

The original x, y and z coordinates can be aligned with the vectors mi, m2 and 
n by first rotating the axes counterclockwise by an angle in the x-y plane (thus 
aligning the new y-axis with m2), then rotating counterclockwise by an angle 6 
about m2. 
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Given that h'^f^ is orthogonal to n, the perturbation tensor can be decomposed 
as the linear combination 



h'.f. = Pmijmik + Q{mijm2k + "^2j"^lA:) + Rm2jm2k, (4.27) 



where 



P = h'^,m{m1 = - ^ (ki - /Co), (4.28) 

vrr sm 6 

Q = h'm{m^2 = ^^(^1 - ^o) (4.29) 



and 



R = h'^.mim', = _ (4.30) 

■' Tcr sm (J 

Thus P = —R and is conventionally denoted as h^, while Q is usually denoted 
as hx- The linearly polarized waveforms therefore have the very simple form^ 

= Jj^h^[(r - t + tL)AL][cos(ri + - 0) - cos(K, + ^ - 0)] (4.31) 

and 

hx = -^^[(r - t + ti) Al] [sin(l^i + ^ - 0) - sin(K) + - </>)], (4.32) 

TXT' sm 1/ 

where it should be recalled that tp = — (k/x)"^ ln[(r — t + tL)/tL\ is a known 
function of the observer time t, while Vi and Vq are implicit functions of t, 6' and 
determined by (I4.1ip and ( 14.121) . 

The expressions (14.311) and fl4.32p for and hx were of course generated on 
the assumption that the observer's y-axis (which is 1112) is parallel to the x-y plane 
of the loop. In the more general situation, the observer's x- and y-axes will be 



^Note that a much quicker and direct way of calculating and hx is to introduce the 4- 
vectors m° — [0, mi]° and 7713 = [0, 1x12]°. Since m'^na = m2na = it follows immediately from 



(|43|) that 

h+ = habmlml 

and 

hx = habm1m2. 
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rotated by some angle ^9 counterclockwise about n from the positions of mi and 
m2. In such the waveform matrix will have the form 



cos {} — sin -(9 
sin -d cos ^ 
1 







-h, 










COS — sin "Q 
sin •& cos d 
1 



/i+ cos 2'i9 — /ix sin2'(9 /ix cos2'(9 + sin2'(9 
/ix cos 2?? + sin 2'(9 — /i^. cos2'(9 + /ix sin2'(9 




and so the linearly polarized waveforms are 



= h+ cos 2d — hx sin 2d 



and 



h'x = hx cos 2d + h+ sin 2'd. 



(4.33) 



(4.34) 



The factor of 2 multiplying the rotation angle d is just a reflection of the fact that 
the graviton is a spin-2 particle. 



5. Plotting the Waveforms 

All that remains now is to plot the linearly polarized waveforms /i+ and /ix as 
functions of T = (t — r)/ti for some sample values of the observer latitude and 
rotation angle d. Note from equations (14.311) and (14.321) that both waveforms are 
bounded by an envelope that falls off linearly with T, but are otherwise periodic 
m. tfj = — (k/x)"^ ln[(r — t + tL)/tL\ with period 27r (as follows from the fact Vq 
and Vi are also periodic in ip). Since different choices of the observer longitude 
simply alter the phase of the waveforms within the envelope, the shapes of the 
waveforms are easily generalized from an examination of the case = only. 

The most problematic part of the calculation is solving the equations (14.111) 
and (I4.12P for Vq and Vi. A natural method of solving equation (14.111) is to use 
the Newton-Raphson algorithm, which involves iterating the recurrence relation 



Jn+l 



Vn — cos(f „ + ■?/') sin 6^ + I cos 6 
1 + smivn + ^j) sin 



(5.1) 



with vq = — |cos6'. If the sequence {vn} converges, it will converge to Vq. Note 
from (14.111) and (14. 12p that if Vq is known as a function of ip then Vi{iIj) = 
—71 — Vo(— V') (on the assumption that cj) = 0). 
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In the particular case Q = -n jl the recurrence relation fl5.ll) does not converge, 
as the denominator is not bounded away from zero. An alternative is to solve for 
Vo by simple iteration, using the recurrence relation 

= cos(t;„ + V^) (5.2) 

with fo = 0. This recurrence relation is of course not a contraction mapping, 
and it seems not to converge without a little tweaking (entering instead a limit 
cycle with period 2). Fortunately, it does seem to converge if fs is set equal 
to \{v2 + "Wi), rather than calculated from (15. 2p . and the algorithm continued as 
before from there. 

For small values of 6*, it is possible to develop power series solutions of equations 
f UrTTD and fHTT^ for Vq and Vi, which when = read: 

\/o = -| + e sin ^/^ + e^{\Ti + cos V sin il)) + 0{f) (5.3) 

and 

V-^ = + e sinip - e^^Ti - costpsin^) + 0(9^). (5.4) 

The corresponding asymptotic expansions for the waveforms and therefore 
describe, to leading order in 6, modulated cosine and sine waves in ip: 

h+ = ^ [{r-t + trMti] cos + 0{e^) (5.5) 

and 



hx=^[{r-t + tL)/tL] sin./; + 0{e^). (5.6) 

Note in particular that the amplitudes of both waveforms tend to zero as 6' — 0, 
as is to be expected from the fact, mentioned earlier, that there is no flux of 
gravitational energy along the z-axis. 

In Figure 3a to 3d the waveforms and /ix (which by definition have a 
rotation angle = 0), and the rotated waveforms h'j^ and h'^ with d = vr/S, are 
plotted in units of ^L^/r for X varying between (on the null surface t = r) and 
1 (on the future light cone of the evaporation point at t = t^,) for an observer 
latitude of 6* = tt/6. Note that the vertical scale is the same on all four plots, so 
that the amplitudes of the waveforms can be compared by eye. 

Figures 4a to 4d are similar, except that now the observer latitude is ^ = 7r/4. 
The vertical scale in Figures 4a to 4d is larger than the vertical scale in Figures 3a 
to 3d by a factor of about 1.33 to 1. Figures 5a to 5d are again similar, except with 
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9 = 71 /3. The vertical scale in Figures 5a to 5d is larger than that in Figures 4a 
to 4d by another factor of about 1.13 to 1. It is evident that the peak amplitude 
of hx increases monotonically as 6 varies from n/6 to 7r/3, but that the peak 
amplitude of /i+ remains roughly constant over this range - a reflection of the fact 
that /i+ — s> at both latitude extremes {9 = and 9 = vr/2), in accordance with 
(EBD and g32D. 

Figure 6b plots for 9 = n/2. Because like cos 6* near the equatorial 

plane, Figure 6a plots not /i+ itself but \img^T,/2{h+/ cos 9). The vertical scale in 
Figures 6a and 6b is about 1.11 times larger than in Figures 5a to 5d, and overall 
about 1.67 times larger than in Figures 3a to 3d. The rotated waveforms h'_^ and 
h'^ are not shown in the case 9 = 71/2 because they are just —^h^ and -^^x, 
respectively. 

Apart from the changes in amplitude, the most obvious trend in the shape 
of the waveforms is the departure from the simple modulated cosine and sine 
waves visible at small values of 9 to the noticeably spiked (in the case of h^) and 
bulbous (in the case of h^/ cos 9) waveforms near 9 = 7i/2. Allen and Ottewill 
[H] have plotted the linearly polarized waveforms from the stationary AGO loop 
over a single cycle for the latitude angles 9 = 7i/2, 7r/4 and 7r/20, and their results 
are recognizably the same as those shown here, although (of course) Allen and 
Ottewill did not include the frequency acceleration induced by the logarithmic 
dependence of ip on t, or the linear modulation factor (r — t + tL)/tL- Oddly, 
Allen and Ottewill's plot of sin0 does not vanish at 9 = 7c/2, even though 
they state explicitly in the text that it should. This discrepancy is presumably 
an artefact of their reliance on truncated Fourier series to plot the waveforms. 

For the sake of visual clarity, the parameter kh appearing in the expression 
for has in Figures 3, 4, 5 and 6 been assigned the value 10^^. This is of course 
not a realistic value for nfi. Given that k ^ 3. 103 7 for the AGO loop, and fi 
is expected to be of order 10~^ or smaller for a cosmic string formed at a GUT 
symmetry-breaking phase transition, the roughly half-dozen cycles visible in each 
of the graphs in Figures 3, 4, 5 and 6 should be replaced by at least 200, 000 
cycles. 

6. Detectability 

The lifetime = -^^;^Lq of an evaporating AGO string loop obviously depends 
on the values assumed for /i (which with dimensional units restored is G/i/c^) 
and the loop's initial length Lq. If they were ever present in the early Universe, 



16 



GUT cosmic strings would have condensed at about 10~^^ seconds after the Big 
Bang. It is also thought that, to a first approximation, the string network will 
have evolved towards a scaling solution, in which long strings are straight on scales 
smaller than a characteristic length comparable to the horizon size, and so loops 
- which form by the intersection or self-intersection of long strings - first appear 
at the horizon scale. 

As an example, AGO string loops with fi ~ 10~^ that formed at the end of 
the radiation-dominated era - about 4000 years after the Big Bang - would have 
had Lq ~ 4000 light years and lifetimes of about 10*^ years. With the waveforms 
executing a total of 200, 000 cycles, this corresponds to a frequency of about 10~^ 
Hz at the time of formation, although subsequent redshifting of the gravitons 
by a factor of (where z^q 1.2 x 10^ for a universe with critical density 
and a Hubble constant of 73 km/s/Mpc |TT]) would reduce this to a currently 
observed frequency of about 10"^'^ Hz. At an observer distance of r ~ 10^° 
light years, the units fiLo/r along the vertical scale in Figures 3, 4, 5 and 6 
would be 4 X 10~^^. Loops formed at earlier epochs and still extant at the end 
of the radiation-dominated era would have been smaller, and their frequencies 
correspondingly higher and astrophysically more interesting (as millisecond pulsar 
timing is sensitive to frequencies of about 10~^ Hz, while the LIGO frequency 
window extends from about 10 to 10^ Hz), but their amplitudes would have been 
proportionally smaller, given that /i+ and h^. scale as L{t). 

However, it is not the waveform amplitudes themselves that determine the 
detectability of a source of gravitational radiation, but rather the density of grav- 
itational energy Pg per logarithmic frequency interval, or spectral density, which 
is conventionally normalized as fig(/) = p~^{fdpg/df), where is the critical 
energy density of the Universe. An upper bound on the magnitude of Qg for a 
single string loop can be estimated as follows. 

A loop forming at a time tp after the Big Bang would have an initial length 
Lq ~ dp. If the loop radiates a power P over a small time interval At cen- 
tered on a time > tp, then at the current time to 3> the total gravita- 
tional energy output PAt would be distributed over a spherical shell with volume 

assuming a spatially-flat Robertson- 
Walker spacetime with scale factor a{t). The corresponding mean density of the 
shell is Pg = PAt/AV, and so if a is assumed to have a power-law dependence 
a{t) cc with any index k e [|, |] an order-of-magnitude estimate of dp^/dt^ 
turns out to be: 

dpJdtRr^c"^Pz^H^%^ (6.1) 
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(on making the replacement a(to)/a(tB) ~ -Zr)- 

Furthermore, if 7° is the radiative efficiency of the loop then its length at time 
t-R is L(tii) ~ ctp — ^^{Gfi/cjitji — tp), and the corresponding frequency at the 
time of emission is /r ~ 2c/L(tR). The currently observed emission frequency 
is therefore / = a{tji)f^/a{to). Since a{tji) cc tj^^^ for emission times during the 
radiation-dominated era, this leads to the estimate 

f-'df/dtn. ~ tR^I + 7°(G'/i/c)tR/L(tR)], (6.2) 

where the term in square brackets is bounded below by | + 7°(G/i/c^) ~ | (this 
bound being achieved when = tp). 
So 

fdpjdf = if-'df/dt^)-\dpjdt^) < 2c-'Pz^' to ^ (6.3) 

where (with units restored) P = ■j'^{Gfj,/c^)'^c'/G. Also, a useful formula relating 
the redshift factor z at time ^r to the current density p^ad of thermal radiation is 

which when substituted into (16. 3p gives: 

n,{f) ^ p:\fdpjdf) < '-f^\Gp/^fn,,^zi{t^/to)\ (6.5) 

where f2rad = Prad/ Pel ^ parameter whose value is currently estimated to be 4.6 x 

10-5 [n]. 

In particular, if ~ 4000 yr - the time of radiation/matter pressure equilib- 
rium - then with z^^ ~ 1.2 x 10^ and to ^ 1.3 x 10^° yr, the value of the spectral 
density fig for an AGO loop with 7° ^ 40 and Gp/c^ ~ 10"^ is bounded above 
by 2 X 10"^^, and scales as p^. By way of comparison, the spectral density of the 
entire stochastic cosmic string background is estimated for Gp/(? ~ 10"^ to have 
the lower bound fig > 1.4 x 10"^, independently of the observed frequency / for a 
wide range of frequencies [13j. The space-based interferometer LISA is projected 
to have a peak sensitivity of fig ~ 10 at a frequency of 10~^ Hz [H], while 
the estimated strain sensitivity of Advanced LIGO (due to go on-line in 2013) at 
frequencies around 100 Hz corresponds to fig ~ 3 x 10~^° [13] 

There is of course nothing peculiar about either the frequencies or amplitudes 
of the waveforms emitted by evaporating AGO loops, save that their lifetimes 
are longer, their waveform amplitudes are typically smaller, and their spectra 
presumably more stable than those of other loop configurations with the same 
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initial length. More detailed and realistic estimates of the effects and detectability 
of gravitational radiation from them can therefore be inferred from similar studies 
in standard references (for example, [12l[ini[l3llllllll])- The current consensus is 
that the strain sensitivity of LIGO is currently one or two orders of magnitude too 
large to detect any stochastic background from cosmic string loops with G/i/c^ > 
10~^, but that millisecond pulsar spindown rates have already ruled out relic 
strings with G/i/c^ > 10~^. It is therefore unlikely that individual string spectra 
will ever be visible to current or foreseeable gravitational wave observatories. 

7. Conclusions 

As has already been mentioned in [3] , the AGO loop is one of the most important 
flat-space cosmic string solutions, as it is possibly the longest-lived of all loop 
configurations, and is moreover the only loop solution (bar the collapsing circular 
loop) whose evolution is known to be analytically tractable. The fact that the 
linearly polarized waveforms and /ix emitted by the evaporating AGO loop 
can be represented by the simple exact forms fl4.3ip and fl4.32p is an added bonus. 

The actual shapes of the waveforms, as plotted in Figures 3, 4, 5 and 6, are 
admittedly not very surprising, as Allen and Ottewill have already calculated 
and plotted the linearly polarized waveforms from the stationary AGO loop in 
[6]. However, their formulas for /i+ and h-^ involved Fourier series rather than 
closed-form expressions and did not include the effects of frequency acceleration 
or amplitude modulation. What has been presented here is therefore a more 
complete treatment of the problem for this particular, and very significant, loop 
configuration. 
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LIST OF FIGURE CAPTIONS: 



Figure 1: y-z projection of the AUen-Casper-Ottewill loop at times r — £ = 0, 
L/16, L/8 and 3L/16 (top row) and r - £ = L/4, 5L/16, 3L/8 and 7L/16 
(bottom row), where the time offset e is 0.02L. The string has been artificially 
thickened for the sake of visibility, and the z-axis is also shown. 

Figure 2: Schematic representation of the evaporation of the AGO loop. The 
thickened line corresponds to the outer envelope of the loop, and Fq and Fl 
are the future light cones of the origin and the evaporation point respectively. 

Figure 3a: The waveform /i+ (in units of /iLq/t) plotted against X for 6* = 7r/6 
and 1? = 0. 

Figure 3b: The waveform (in units of /iLo/r) plotted against T for ^ = 7r/6 
and = 0. 

Figure 3c: The waveform h'^ (in units of /iLo/r) plotted against X for ^ = 7r/6 
and Ti) — 7r/8. 

Figure 3d: The waveform /i'^ (in units of /jlLo/v) plotted against X for ^ = 7r/6 
and 1? = 7r/8. 

Figure 4a: The waveform h+ (in units of iiLq/t) plotted against T for 6* = 7r/4 
and 1? = 0. 

Figure 4b: The waveform /ix (in units of /iLq/t) plotted against X for 6* = 7r/4 
and 1? = 0. 

Figure 4c: The waveform h'_^ (in units of /iLo/r) plotted against X for ^ = 7r/4 
and "d — Tr/8. 

Figure 4d: The waveform h'^ (in units of /iLo/r) plotted against X for ^ = 7r/4 
and i9 = tt/S. 

Figure 5a: The waveform h+ (in units of iiLq/t) plotted against X for 6' = 7r/3 
and ^? = 0. 

Figure 5b: The waveform (in units of iiLq/t) plotted against X for 6* = 7r/3 
and 1? = 0. 
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Figure 5c: The waveform h'j^ (in units of fJ^Lg/r) plotted against T for ^ = 7r/3 
and 1? = 7r/8. 

Figure 5d: The waveform h'^ (in units of iiLq/t) plotted against T for ^ = 7r/3 
and ^ = 7r/8. 

Figure 6a: The rescaled waveform hj^/ cos 9 (in units of iiLq/t) in the limit as 
^ — >• 7r/2 plotted against T for = 0. 

Figure 6b: The waveform hy, (in units of iiLo/r) plotted against % ior 9 — 7r/2 
and "d — 0. 
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Figure 3b 
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